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ABSTRACT 

We present accurate positions for 857 sources derived from the astrometric analysis of 16 eleven-hour ex- 
periments from the Very Long Baseline Array imaging and polarimetry survey at 5 GHz (VIPS). Among the 
observed sources, positions of 430 objects were not previously determined at milliarcsecond-level accuracy. 
For 95% of the sources the uncertainty of their positions ranges from 0.3 to 0.9 mas, with a median value of 
0.5 mas. This estimate of accuracy is substantiated by the comparison of positions of 386 sources that were 
previously observed in astrometric programs simultaneously at 2.3/8.6 GHz. Surprisingly, the ionosphere con- 
tribution to group delay was adequately modeled with the use of the total electron content maps derived from 
GPS observations and only marginally affected estimates of source coordinates. 
Subject headings: astrometry — catalogues — surveys 



1. INTRODUCTION 

At high galactic latitudes 70% (662 of 1043) of the 7- 
ray bright sources detected in the LAT first-year catalogue 
(1FGL; Abdoetal. (2010)) are associated with AGN with 
high confidence (P > 80%). The vast majority of the Fermi 7- 
ray sources are blazars, with strong, compact radio emission. 
These blazars exhibit fiat radio spectra, rapid variability, com- 
pact cores with one-sided parsec-scale jets, and superluminal 
motion in the jets (Marscher et al. 2006, Taylor et al. 2007, 
Linfordet al. 2011). In the zone \b\ > 10° radio emission at 
parsec scales has been detected with VLBI for 60% (624 out 
of 1040) of the 1FGL 7-ray sources (Kovalev & Petrov, 2011 
in preparation). 

In preparation for the Fermi mission, an all-sky catalogue of 
compact radio sources, the Combined Radio All-sky Targeted 
Eight gigaHertz Survey (CRATES; Healey et al. (2007)) con- 
taining over 11,000 sources was established. This cata- 
logue is being used with considerable success in associat- 
ing Fermi detections with AGN (Abdoet al. 2010). A sub- 
set of 1127 sources from CRATES was observed with VLBI 
to form the VLBI Imaging and Polarimetry Survey (VIPS; 
Helmboldt et al. (2007)). 

The scientific goals for VIPS are: 

1. To gather statistically complete information on large 
subsamples in order to study the apparent properties of 
objects as a function of the angle between the jet and 
the line-of-sight. 

2. To explore the underlying patterns in the jets, which can 
be obscured in images by the peculiarities of individual 
objects. In a manner similar to weak lensing, it is pos- 
sible to average out these peculiarities. 

3. To explore the evolution of AGN by studying the size 
(and eventually age) distribution of Compact Symmet- 
ric Objects (CSOs); (Tremblay et al.(201 1), in prepara- 
tion). 

Leonid.Petrov@lpetrov.net 
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4. To search for additional Supermassive Binary Black 
Hole systems (SBBHs). These can be used to probe 
the stalling radius for black hole aggregation, and to 
determine the prevalence of mergers that might be 
sources of gravitational radiation. The same observa- 
tions will also be used to detect or improve the limits 
on milliarcsecond-scale gravitational lenses. 

The original analysis of VIPS observations was made us- 
ing an automated approach developed by Taylor (2005) based 
on the software package Difmap. That pipeline produces im- 
ages in Stokes parameters /, Q, and U that were consecutively 
analyzed for deriving total intensity, polarized intensity and 
electric vector position angle. These results are discussed 
in depth in Helmboldt et al. (2007). Originally, absolute as- 
trometry was not considered as a goal of this project. Al- 
though it was feasible to reduce VLBA data with AIPS and 
determine source coordinates in the absolute astrometry mode 
(Petrov et al. 2009), that approach raised the detection limit 
by a factor of -\/8 = 2.83, required intensive manual work, 3-5 
working days per experiment, and was not considered practi- 
cal for large projects. 

During 2005-2011, significant efforts were made to de- 
velop modern algorithms for processing survey-style experi- 
ments and to implement them in software. This approach was 
oriented on unattended astrometric analysis of an arbitrarily 
large amount of data and has been described in complete de- 
tail in Petrov et al. (2011). The validation of that approach 
against a set of 500,000 VLBA observations under regular 
geodesy and absolute astrometry program RDV previously 
processed with AIPS has proven its vitality, high accuracy, 
lack of systematic errors, and a low detection limit approach- 
ing the theoretical limit. The amount of time required for pro- 
cessing a typical survey experiment was reduced to 2-3 hours, 
making it practical to re-analyze old VLBA experiments and 
to derive positions of observed sources. 

In this paper we show the results of our re-analysis of VIPS 
observations aimed at absolute astrometry and present the cat- 
alogue of source positions. The motivation of this work was 
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first to add milliarcsecond-level accurate position information 
to each data-sheet associated with the VIPS sample: images at 
5 GHz, polarization measurements, optical photometry from 
the Sloan Digital Sky Survey, redshifts, and 7-ray fluxes. Our 
second goal was to extend the VLBA Calibrator list. Since 
it was known from image analysis that the majority of VIPS 
sources are compact, we expect they will be excellent phase- 
calibrators as well. Our third goal was to obtain highly ac- 
curate positions of optically bright quasars from the VIPS list 
for validation of results from the planned space astrometry 
mission (Lindegren et al. 2008). 

In section 2 we briefly describe the observations, sample 
definition, and scheduling. In section 3 we outline the data 
analysis methodology. The catalogue is presented in section 
4 followed by a brief discussion. Final remarks are given in 
the summary. 

2. OBSERVATIONS 

2.1. Sample Definition 

To meet the primary goals of VIPS, a relatively large sam- 
ple of likely AGN, preferably with data from other wavelength 
regimes, was required. To this end, the Cosmic Lens All-Sky 
Survey (CLASS; Myers et al. 2003) as the parent sample was 
chosen. CLASS is a VLA survey of ~12,100 flat-spectrum 
objects (a > -0.5 between 4.85 GHz and a lower frequency), 
making it an ideal source of likely AGN targets to be fol- 
lowed up with the VLBA. The sample was also restricted to 
lie within the survey area, or "footprint" of the Sloan Dig- 
ital Sky Survey (SDSS; Yorket al. 2000) as defined at the 
outset of VIPS in the fifth data release of the SDSS (DR5; 
Adelman-McCarthy et al. 2006). In DR5 the imaging covers 
8,000 square degrees and includes ~2x 10 s objects. Spec- 
troscopy was obtained as part of the SDSS for ~ 10 6 of these 
objects, about 10 5 of which are quasars. The source cata- 
log was chosen such that all sources lie on the original SDSS 
footprint with an upper declination limit of 65° (imposed to 
avoid the regions not imaged through DR5). We also excluded 
sources below a declination of 15° because it is difficult to 
obtain good (m, v) coverage with the VLBA for these objects. 
To keep the sample size large but manageable and to obtain 
a high detection rate without phase referencing, all CLASS 
sources were selected within this area on the sky with flux 
densities at 8.5 GHz greater than 85 mJy, yielding a sample of 
1,127 sources. For further details see Helmboldt et al. (2007). 

2.2. VIPS Observations 

The VIPS program was allocated 200 hours of observing 
time with the VLBA. This time was divided into 18 observ- 
ing sessions of ~1 1 hours each which ran between January 3 
and August 12 2006 (project BT085). In each session a group 
of 52-54 target sources were observed along with the calibra- 
tors 3C279, and J 1 3 1 0+3220, and two other calibrators drawn 
from the following list: DA193, OQ208, 3C273, J0854+2006, 
and Jl 159+2914. In total 958 VIPS sources were observed. 
Each VIPS target was observed for approximately 500 sec- 
onds divided into 10 separate scans. All observations were 
conducted with four 8 MHz wide, full polarization frequen- 
cies centered at 4609, 4679, 4994, and 5095 MHz. These fre- 
quencies were chosen in order to allow for determination of 
Faraday rotation measures in those sources with sufficiently 
polarized components, and in order to improve the (w, v) cov- 
erage obtained. 

All VLBA observations were scheduled using version 6.05 



of the SCHED program 1 . Using built-in data regarding the 
locations and operation of the VLBA stations, a special op- 
timization mode "HAS" of SCHED automatically produced 
a schedule for a list of targets with scan durations, a starting 
local sidereal time (LST), and total experiment duration that 
is optimized both for (u,v) coverage and efficiency. For each 
observing run, the starting LST and scan time per source was 
varied to produce a schedule that obtained for the entire dura- 
tion of 1 1 hours the vast majority (if not all) of the required 
scans for a subset of ^50 VIPS targets while minimized slew- 
ing time between sources. 

Scan lengths were set to ~50 seconds, and the scheduling 
software distributed 10 scans of each target source at approxi- 
mately equal hour angle intervals. Optimization was done for 
the VLBA antenna PIETOWN since it is fairly centrally located 
to the array. The minimum elevation at that antenna was set 
to 15 degrees, and the minimum number of antennas above 
the horizon was set to 8. The maximum deviation from the 
optimal time for the scan was set at 3 hours. 

Care was also taken to select the correct polarized calibra- 
tion source(s) for each run so that it/they would be observed 
over a wide range of parallactic angle values while not signif- 
icantly reducing the efficiency of the schedule for that run. 
These sources were also used as amplitude calibrators for 
evaluation of antenna complex bandpasses. 

The schedule optimization that minimized slewing time re- 
sulted in scheduling sources at a relatively wide band across 
azimuths and elevations at observing stations. See Figure 1 as 
an example. No elevation cutoff-limit was enforced, and an- 
tennas were used as low as the physical horizon mask permits, 
down to 3°. 
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Figure 1. Azimuthal-elevation source distribution over the sky at station HN- 
VLBA during 11-h long experiment bt()85f. Source 1609+ 190 is highlighted 
by a big green rilled circle. Other sources are displayed as small outlined grey 
circles. 

Unfortunately, the archive magnetic tape with correlator 
output from two experiments, bt085d and bt085e, had dete- 
riorated to the extent that it could not be read, and the data 

1 See documentation at 
http : //www . aoc . nrao . edu/sof tware/sched 
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were completely lost. Therefore, we were able to re-analyze 
only 16 out of 18 experiments for a total of 857 sources. 

3. DATA ANALYSIS 

The data were correlated at the Socorro hardware VLBA 
correlator. The correlator computed the spectrum of cross cor- 
relation and autocorrelation functions with a frequency reso- 
lution of 0.5 MHz and accumulation intervals 1.96608 s long. 

3.1. Post-correlation analysis 

The procedure of further analysis is described in complete 
detail in Petrov et al. (2011). Here only a brief outline is 
given. First, the fringe amplitudes were corrected for the sig- 
nal distortion in the sampler and then were calibrated accord- 
ing to measurements of system temperature and elevation- 
dependent gain. Phase of the phase-calibration signal was 
subtracted from the fringe phases. Then the group delay, 
phase delay rate, group delay rate, and fringe phase were de- 
termined for the RR polarization (i.e. right circular polariza- 
tion (RCP) at a reference station correlated against the RCP 
at a remote station of a baseline) for all observations, us- 
ing the wide-band fringe fitting procedure. These estimates 
maximize the sum of the cross-correlation spectrum coher- 
ently averaged over all accumulation periods of a scan and 
over all frequency channels in all intermediate frequencies 
(IFs). After the first run of fringe fitting, an observation at 
each baseline with the reference station (KP-VLBA) with the 
strongest signal to noise ratio (SNR) were used to compute the 
station-based complex bandpass corrections of the RR cross- 
correlation spectrum. The initial bandpass is the cross spec- 
trum coherently averaged over time divided by the average of 
each IF. If the bandpass characteristics of each IF were con- 
stant over an experiment, then applying the bandpass, i.e. di- 
viding the cross-spectrum by the bandpass, would make the 
cross-spectrum averaged over time completely fiat, which re- 
duces the coherence loss. This procedure was repeated for 
12 sources with the strongest SNR at each baseline, and the 
results were averaged. 

Next, we computed the polarization complex bandpass. 
This was done for each baseline with a reference station in 
two steps. First, we selected a scan with the highest SNR, 
applied the initial RR complex bandpass and computed group 
delay r g , group delay rate f g , phase delay t p , and phase delay 
rate f p at the reference frequency using RR data. Then we ap- 
plied the RR complex bandpass to the LL data CS at the kth 
accumulation period at the y'th spectral channel and rotated its 
phase using parameters of the fringe fitting to the RR data: 

QLL x e -j'(oio(Tp + f p (t k -t )) - 2(l/>i-^ 2 ) + (,Wj-WQ)T g ) ^ 

where wo is the reference frequency, to is the fringe reference 
time, ojj is the frequency of the y'th spectral channel, and % is 
time of the kth accumulation period. We also subtracted the 
doubled difference of parallactic angles (ipi — ip2) at the refer- 
ence and remote stations of a baseline. The latter correction 
accounts for the fact that the phase rotation due to a change 
of the receiver orientation with respect to some reference par- 
allactic angle is opposite for the RCP and LCP (left circular 
polarization) signal. The residual amplitudes were normal- 
ized at each IF by its average over all spectral channels. We 
approximated the normalized residual complex spectrum with 
the Legendre polynomial of the 5th degree. These polynomi- 
als define the so-called initial polarization bandpass. It was 



found that using polynomials of higher degree did not provide 
further improvements. 

Next, we processed 12 observations with the highest SNR 
at each baseline with the reference station. We applied the 
sampler correction, phase calibration to the RR data, com- 
puted parameters of the fringe fitting to the RR data, applied 
the sampler correction, phase calibration, performed the phase 
rotation according to expression 1 to the LL data of the same 
observation, and divided the residual spectrum by the initial 
complex polarization bandpass. The phase and the logarithm 
of the amplitude of the residual LL spectrum are used as 
the right-hand side of the observation equations that approx- 
imate the residual complex polarization bandpass with Leg- 
endre polynomials of the 5th degree polynomial for each IF 
and station, except the reference station. This system of lin- 
ear equations is solved using least squares (LSQ). The product 
of these estimates with the initial polarization bandpass forms 
the final complex polarization bandpass. 

The procedure of computing group delays was repeated 
with the complex polarization bandpass applied to all the data 
and using the Stokes parameter /, i.e. the linear combina- 
tion of RR and LL cross-spectrum Crr and Cll- 1/2(C«« + 
Cll/P^/Brr, where Brr is the complex bandpass for the RR 
data, and Pb is the polarization bandpass. 

If the electronics of every IF were stable during the ex- 
periment and observed sources are completely unpolarized, 
the SNR of the coherent sum over time and frequency of the 
Stokes parameter / would have been a factor of V2 times the 
SNR when only RR or LL polarization data are used. Anal- 
ysis of the SNR ratios of Stokes parameter / data to the SNR 
of RR data showed that achieved SNR was on average 96% of 
its theoretical value. 

This part of the analysis is done with V1A4A software 2 . 

3.2. Astrometric analysis 

The results of fringe fitting the linear combination of RR 
and LL data were exported to the VTD/post-solve VLBI anal- 
ysis software 3 for interactive processing of the group delays 
with SNR high enough to ensure the probability of false detec- 
tion is less than 0.001. This SNR threshold is 5.5 for the VIPS 
experiment. A detailed description of the method for evalu- 
ation of the detection threshold can be found in Petrov et al. 
(2011). Next, theoretical path delays were computed accord- 
ing to the state-of-the art parametric model as well as their 
partial derivatives. Small differences between group delays 
and theoretical path delay were estimated using a least squares 
adjustment to the parametric model that describes the obser- 
vations. Coordinates of target sources, positions of all sta- 
tions (except the reference station), parameters of the spline 
that models residual path delay in the neutral atmosphere in 
the zenith direction for all stations, and parameters of another 
spline that models the clock function, were solved for in LSQ 
solutions that used group delays of the Stokes parameter /. 
Outliers of the preliminary solution were identified and tem- 
porarily barred from the solution. 

The most common reason for an observation to be marked 
as an outlier is a misidentification of the main maximum of 
the two-dimensional Fourier-transform of the cross-spectrum. 
The procedure of fringe fitting was repeated for observations 
marked as outliers. Using parameters of the VLBI model ad- 
justed in the preliminary LSQ solution, we can predict group 

2 Available at http : //astrogeo . org/pima. 

3 Available at http : //astrogeo . org/vtd. 
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delay for outliers with accuracy better than 1 ns. Then we 
re-run the fringe fitting procedure for observations marked as 
outliers in the narrow search window ±1 ns within the pre- 
dicted value of the group delay. New estimates of group de- 
lays with probabilities of false detection less than 0.1, which 
corresponds to the SNR > 4.6 for the case of a narrow fringe 
search window, are used in the next step of the preliminary 
analysis procedure. The observations marked as outliers in the 
preliminary solution and detected in the narrow window at the 
second round of the fringe fitting procedure, were tried again. 
If the new estimate of the residual was within 3.5 standard 
deviations, the observation was restored and used in further 
analysis. Parameter estimation and elimination of remaining 
outliers were then repeated. Finally, the additive weight cor- 
rections were computed in such a way that these corrections, 
added in quadrature to the a priori weights computed on the 
basis of group delay uncertainties evaluated by the fringe fit- 
ting procedure, will make the ratio of the weighted sum of 
residuals to their mathematical expectation close to unity. 

In total, 310653 estimates of group delays out of 347803 
scheduled were used in the solution. The number of observa- 
tions of 95% of the target sources was in the range 200-440, 
with the median number 377. Each source was detected, ex- 
cept 1058+245. Detailed analysis showed that there are two 
objects in the field of view near 0748+582. The companion 
0748+58A is 72". 8354 apart. 

The result of the preliminary solution provided the clean 
dataset of group delays with updated weights. We used this 
dataset and all dual-band S and X data acquired under abso- 
lute astrometry and space geodesy programs from April 1980 
through April 201 1, in total 8.2 million observations, to obtain 
final parameter estimation. Thus, the VIPS experiments were 
analyzed by exactly the same way as all other VLBI exper- 
iments. The estimated parameters were right ascensions and 
declination of all sources, coordinates and velocities of all sta- 
tions, coefficients of the B-spline expansion of non-linear mo- 
tion for 18 stations, coefficients of harmonic site position vari- 
ations of 48 stations at four frequencies: annual, semi-annual, 
diurnal, semi-diurnal, and axis offsets for 69 stations. Esti- 
mated variables also included the Earth orientation parameters 
for each observing session and parameters of clock function 
and residual atmosphere path delays in the zenith direction 
modeled with the linear B-spline with interval 60 minutes. 
All parameters were adjusted in a single LSQ run. 

The system of LSQ equations has an incomplete rank and 
defines a family of solutions. In order to pick a specific el- 
ement from this family, we applied the no-net rotation con- 
straints on the positions of 212 sources designated as "defin- 
ing" in the ICRF catalogue (Ma et al. 1998) that required the 
positions of these sources in the new catalogue to have no ro- 
tation with respect to the position in the ICRF catalogue. No- 
net rotation and no-net translation constraints on site positions 
and linear velocities were also applied. The specific choice of 
identifying constraints was made to preserve the continuity of 
the new catalogue with other VLBI solutions made during last 
15 years. 

The global solution sets the orientation of the array with 
respect to an ensemble of ~5000 extragalactic remote radio 
sources. The orientation is defined by the continuous series 
of Earth orientation parameters and parameters of the empir- 
ical model of site position variations over 30 years evaluated 
together with source coordinates. Over 400 common sources 
observed in the VIPS provided a firm connection between the 
new catalogue and the old catalogue of compact sources. 



3.3. Accounting for the ionosphere contribution 

Since the observations were made at a relatively low fre- 
quency of 5 GHz, one may expect that the dominant source of 
errors will be the mismodeled ionosphere contribution. The 
superior way to alleviate the frequency-dependent contribu- 
tion of the ionosphere to group delays is to record simulta- 
neously at two widely separated frequencies, i.e. 2.3 GHz 
(S-band) and 8.6 GHz (X-band) and to use ionosphere-free 
linear combinations of these observables in the data analysis. 
The residual delay due to high order terms that are not elimi- 
nated in the combination of observables is at a level of several 
ps (Hawarey et al. 2005), which is negligible. 

In the absence of simultaneous dual-band observations, we 
computed the ionospheric path delays from the TEC maps de- 
rived from the analysis of Global Positioning System (GPS) 
observations using the analysis center CODE (Schaer 1998). 
Using GPS-derived TEC maps for data reduction of astro- 
nomic observations, first suggested by Ros (2000), has be- 
come a traditional approach for data processing. However, 
we should be aware that these path delays are an approxima- 
tion of the ionospheric contribution and they cannot account 
completely for the contribution of the ionosphere. First, the 
TEC model has rather a coarse resolution: 5°in spatial coor- 
dinates and 2 in time. Therefore, ionospheric fluctuations on 
time scales less than several hours cannot be represented by 
this model. Second, the model itself is not precise, since the 
number of GPS stations that were used for deriving the model 
(179 in January 2006) is relatively small and their distribution 
around the globe is far from homogeneous. 

The availability of dual-band observations at the VLBA net- 
work during the period of time overlapping with the VIPS 
campaign prompted us to use this dataset for characterizing 
the residual errors of ionospheric path delays from GPS TEC 
maps. Since the ionospheric electron density varies by more 
than one order magnitude during a Solar activity cycle, we 
restricted our analysis to 15 twenty-four hour experiments 
under geodetic program RDV (Petrov et al. 2009) and VCS6 
(Petrov et al. 2008) over the period of time from June 2005 
through March 2007. In total, over 130,000 dual-band obser- 
vations were used for comparison. For each observation, we 
computed the contribution of the ionosphere to group delay at 
4.8 GHz from both X-band and S-band VLBI group delay ob- 
servables and from GPS TEC maps. The technique of further 
analysis is described in details in Petrov et al. (2011). Here 
we briefly outline it. 

For each baseline and each experiment, we computed two 
statistics: RV, which stands for root mean square (rms) of the 
differences in the ionospheric group delay between VLBI and 

GPS measurements, and RG = ^Jvms 2 gl +rms^ 2 , where rms,,, 

is the rms of the GPS path delay at the i th station of a base- 
line during a session. As we found previously, these statistics 
are highly correlated. For each baseline we computed a linear 
regression between them: RV/, = Fi,+Sb RG/,. The coefficients 
Fb and Sb depend on the baseline. Using this regression, we 
can predict for a given baseline the rms of the residual contri- 
bution of group delay computed from GPS TEC maps. 

3.4. Error analysis 

One half of the sources scheduled for VIPS has 
been observed as part of the VLBA Calibrator Survey 
(Beasley et al. 2002, Fomalont et al. 2003, Petrov et al. 2005, 
2006, Kovalev et al. 2007, Petrov et al. 2008) and RDV pro- 
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Figure 2. Differences in positions of common sources derived from analysis of VIPS observations at C-band against positions derived from analysis of dual-band 
X/S observations. 



gram (Petrov et al. 2009) aiming to determination of source 
positions in the absolute astrometry mode with the highest 
accuracy. These observations were made simultaneously at 
X and S-bands. We can compare the positions of common 
sources from single-band VIPS observations with dual-band 
observations and, considering the positions from dual-band 
observations as the ground truth, derive the error model. 

We scrutinized the list of 430 common sources and dis- 
carded from the list 54 objects with position uncertainties 
from X/S observations exceeding 1 mas. We are aware that 
the positions of 386 common sources from the VIPS cata- 
logue were made from both VIPS observations and all others 
observations. We cannot exclude all common sources from 
the VIPS solution, since in that case the matrix of observa- 
tions will be altered and properties of such a solution would 
be significantly different from the properties of the VIPS cat- 
alogue. To circumvent this problem, we sorted the list of re- 
maining 386 sources in ascending order of right ascensions 
and split it into two subsets of 193 sources each, even and odd. 
Then we ran two special solutions: solution A that had 193 
sources from the even subset excluded in all experiments, ex- 
cept VIPS, and solution B that had 193 sources from the odd 
subset excluded in all experiments, except VIPS. The solution 
setup, except the differences in the source list, was identical 
to that used for deriving the VIPS catalogue. In all solutions 
we estimated positions of 5317 sources. We obtained posi- 
tions of common sources from these two solutions derived 
only from the C-band VIPS observations and then combined 
them to form a catalogue of 386 objects. 

A comparison of these two catalogues is shown in Figure 2. 
We present the source position differences in right ascensions 
and declinations as a function of source declinations. We 
see that for a majority of sources, the differences look like 
a random noise and do not exhibit any systematic pattern. 
There are several outliers, the most notable is 0820+225 (See 
Figure 3 and Table 1). Analysis of images 4 of 0820+225, 
1108+201, 1323+321, and eight others revealed a common 
pattern: all are sources with multiple distinctive components, 
and the relative strength of components varies between bands. 
Therefore, source positions are related to different parts of the 
source. A similar phenomenon was noticed during analysis of 
differences of X/S and K-band (24 GHz) positions reported 

4 Brightness distributions in FITS-format as well as images 
for these and all VIPS sources are available on the Web from 
http : //www . phys . unm . edu/~gbtaylor/VIPS/ and from 
http : //astrogeo . org/vlbi_images 



by Petrov et al. (2011). On average, 1% sources exhibit this 
behavior. 
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Figure 3. Image of J0823+2223 (0820+225) in four frequency bands. The 
center of the image is set to the brightest component. The images at X and 
S bands are produced from analysis of VCS 1 observations (Kovalev, private 
communication), the image at C-band was produced from analysis of VIPS 
observations, and the image at U-band was produced from analysis of obser- 
vations from the MOJAVE program (Lister et al. 2009). 



After removing the outliers, we sought for the ad hoc vari- 
ance that, when added in quadrature with the formal error, 
makes the ratio of the weighed sum of the differences to their 
mathematical expectation close to unity. We found that the 
additional variance for right ascensions scaled by cosine of 
declination is 0.24 mas and the variance for declinations is 
0.31 mas. The weighted root mean squares (wrms) of the dif- 
ferences in right ascensions scaled by cosine of declination 
is 0.44 mas and the wrms of the differences in declination is 
0.58 mas. 

We should note that our estimate of the additive variances 
provides an upper limit on the errors. First, de-selecting 
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Table 1 

The table of sources with the greatest differences between VIPS and X/S 
positions. Quoted flux densities are the median correlated flux densities at 
baselines longer than 5000 km. 



200 



Dist 
(mas) 


Err 
(mas 


B1950-name 


J2000-name 


Comment 




14.6 


0.5 


0820+225 


J0823+2223 


Long, complex jet 




5.0 


0.4 


1323+321 


J 1326+3 154 


CSO 




3.6 


0.5 


1721+343 


J1723+3417 


Long, complex jet 




3.4 


0.4 


1048+347 


J 1050+3430 


Double source (core-jet) 


3.4 


0.4 


1108+201 


Jllll+1955 


CSO 




3.1 


0.4 


1109+353 


Jl 112+3503 


Weak: 76 mjy at 8 


GHz 


3.0 


0.4 


1600+244 


J1602+2418 


CSO candidate 




2.9 


0.4 


1626+229 


J 1628+2247 


Long, complex jet 




2.8 


0.4 


1019+309 


J1022+3041 


Long, complex jet 




2.6 


0.5 


1117+543 


Jl 120+5404 


Weak: 38 mjy at 8 


GHz 


2.6 


0.7 


0752+639 


J0756+6347 


Strong core-jet 




2.4 


0.4 


0932+367 


J0935+3633 


Double source (core-jet) 


2.2 


0.5 


1221+464 


J1223+4611 


Weak: 67 mjy at 8 


GHz 


2.1 


0.4 


1525+314 


J1527+3115 


Long, complex jet 





0) 

o 

f-H 

o 

tsi 



sources in trial solutions A and B weakened these solutions 
with respect to the solution used for deriving the VIPS cata- 
logue and increased position errors of common sources. Sec- 
ond, the additive variances 0.24 and 0.31 mas are comparable 
with errors of the reference X/S solution. Therefore, the ad- 
ditive variances reflect contributions of errors of both C-band 
VIPS and the X/S catalogues. 

4. SOURCE POSITION CATALOGUE 

Table 2 displays 12 out of 857 rows of the VIPS catalogue 
of source positions. The full table is available in the elec- 
tronic attachment. Column 1 contains the source flag: "X" 
if the sources was previously observed at X and S-bands un- 
der absolute astrometry programs, "C" otherwise. Columns 
2 and 3 contain the J2000 and B1950 IAU names; column 
4, and 5 contain right ascensions and declinations. Columns 6 
and 7 contain positional uncertainties in right ascension (with- 
out multiplier cos S) and in declinations with noise added in 
quadrature to their formal errors. Column 8 contains correla- 
tion coefficients between right ascension and declination, col- 
umn 9 contains the total number of observations used in posi- 
tion data analysis. Columns 10 through 15 contain the median 
estimates of the correlated flux density over all experiments 
of a given source in two ranges of baseline projection lengths: 
less than 900 km, and longer than 5000 km. The first estimate 
is close to the total correlated flux density, the second estimate 
is close to the correlated flux density of the unresolved com- 
ponent. The estimates of the median correlated flux densities 
are given for three bands: S-band (columns 10,11), C-band 
(columns 12,13), and X-band (columns 14,15). The estimate 
of the flux density at C-band was computed from analysis of 
VIPS data by Helmboldt et al. (2007), estimates at S and X 
band are from re-analysis of VCS and RDV observations. If 
there were no detections at the range of the baseline projection 
lengths, then -1 .000 is put in the table cell. 

The histogram of the semi-major axes of position uncertain- 
ties is given in Figure 4. 

5. DISCUSSION 

The astrometric accuracy of VIPS exceeded our expecta- 
tions. Although absolute astrometry was not considered at 
all during design of the experiment, position uncertainties of 
target sources observed in single band VIPS experiments at 
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Position uncertaintiy (mas) 

Figure 4. The histogram of semi-major axes of inflated position uncertain- 
ties of the VIPS source position catalogue. The last bin shows uncertainties 
exceeding 1 mas. 

5 GHz turned out better than the position uncertainties of 
dedicated VCS experiments. The contribution of the residual 
ionosphere path delay after applying the GPS-derived TEC 
model was small and did not affect the positional accuracy 
significantly. 

Several factors contributed to this success. First, the wide- 
band frequency setup spanned over 494 MHz was used. This 
is essential. If a contiguous bandwidth of 32 MHz had been 
recorded, as is often selected in imaging experiments, group 
delays would have been a factor of 22 less precise, and as 
a result, the source position accuracy would have been more 
than one order of magnitude worse. Second, each source was 
scheduled in 10 scans equally distributed over the hour an- 
gle of the array center (station PIETOWN). The schedule was 
optimized to get good (m, v)-coverage with no elevation cut- 
off enforced. Observations at elevations as low as 3° were 
scheduled. Good (m, v)-coverage also provides a good sepa- 
ration among variables of interest such as source coordinates, 
as well as among nuisance parameters such as the clock func- 
tion. Frequent scheduling of low-elevation observations im- 
proves the reliability of estimates of residual atmosphere path 
delay in zenith parameters, since it helps to decorrelate these 
parameters with station positions and clock function param- 
eters. Third, the integration time was correctly predicted, so 
the vast majority of sources were detected at both short and 
long baselines. A lack of observations at long baselines re- 
duces the accuracy of source positions significantly. Finally, 
the experiments were carried out during the Solar minimum 
when ionosphere disturbances were exceptionally low. 

Of the 14 sources with significant discrepancies between 
their position measured in VIPS and that derived from dual- 
band X/S observations (see Table 1, the majority (10/14) have 
complex structure in which the brightest component is a func- 
tion of angular resolution and/or observing frequency. This 
includes three Compact Symmetric Objects which generally 
have more complex structure than core-jets. Of the remainder, 
three sources are core-jets, but quite weak (< 100 mJy/beam), 
so it is likely that in this case the errors are less systematic 
and more the result of low SNR. The single remaining source, 
VIPS J0756+6347, is a strong and ordinary core-jet and the 
reason for the discrepancy in positions is unclear. 

6. CONCLUSIONS 

We have determined positions of 857 sources observed in 
16 VIPS experiments. Milliarcsecond-accurate coordinates of 



Precise absolute astrometry from VIPS 



7 



Table 2 

The first 12 rows of the VIPS catalogue of source positions of 857 sources. The table columns are explained in the text. The full table is available in the electronic 

attachment. 



Correlated flux density (in Jy) 



Flag 


Source Names 


J2000 Coordinates 


Errors (mas) 




S-band 


C-band 


X-band 




IVS 


IAU 


Right ascension 


Declination 


Aq 


A<5 


Corr 


# Obs 


Total 


Unres 


Total 


Unres 


Total 


Unres 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


X 


0552+398 


J0555+3948 


05 55 30.805615 


+39 48 49.16493 


0.40 


0.24 


-0.032 


1775 


3.457 


2.765 


6.957 


4.226 


4.736 


2.144 


C 


07164450 


J07 19+4459 


07 19 55.511712 


+44 59 06.84827 


1.68 


1.52 


-0.059 


98 






0.151 


<0.02 






c 


0722+393 


J0726+3912 


07 26 04.737260 


+39 12 23.32117 


0.55 


0.57 


-0.300 


302 






0.076 


0.045 






c 


0722+415 


J0726+4124 


07 26 22.420172 


+41 24 43.66844 


0.53 


043 


-0.217 


362 






0.106 


0.090 






c 


0722+615 


J0726+6125 


07 26 51.681732 


+61 25 13.68060 


0.90 


0.42 


-0.017 


336 






0.096 


0.088 






X 


0723+488 


J0727+4844 


07 27 03.100595 


+48 44 10.12669 


0.52 


0.34 


-0.004 


362 


0.447 


0.357 


0.263 


0.222 


0.330 


0.280 


X 


0729+562 


J0733+5605 


07 33 28.615274 


+56 05 41.73927 


2.34 


1.98 


-0.238 


86 


0.099 


1.269 


0.113 


0.022 






c 


0731+595 


J0735+5925 


07 35 56.300840 


+59 25 22.11156 


0.94 


0.48 


0.022 


322 






0.053 


0.037 






X 


0733+300 


J0736+2954 


07 36 13.661088 


+29 54 22.18498 


0.41 


042 


-0.038 


344 


0.464 


0.409 


0.273 


0.158 


0.302 


0.178 


c 


0733+287 


J0736+2840 


07 36 31.198545 


+28 40 36.80515 


2.05 


1.97 


-0.144 


58 






0.039 


0.025 






X 


0733+261 


J0736+2604 


07 36 58.073689 


+26 04 49.94573 


0.41 


0.42 


-0.034 


351 


0.336 


0.086 


0.240 


0.108 


0.326 


0.215 


c 


0734+269 


J0737+2651 


07 37 54.975269 


+26 51 47.44790 


0.54 


0.60 


-0.141 


302 






0.071 


0.051 







Note. — Table 2 is presented in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form and contents. 
Units of right ascension are hours, minutes and seconds. Units of declination are degrees, minutes and seconds. 



430 objects were determined with VLBI for the first time. 

We found that during the Solar minimum the residual iono- 
sphere contribution to the group delay after applying the re- 
duction for path delay from GPS TEC map is small enough 
not to cause systematic errors exceeding 0.3 mas at frequen- 
cies as low as 5 GHz. 

Our error analysis showed that for 95% of the sources the 
uncertainty of their positions range from 0.3 to 0.9 mas. This 
estimate is substantiated by comparison with a large set of 
386 VIPS sources that were also observed in X/S absolute as- 
trometry programs. The position uncertainty of the VIPS cat- 
alogue is better than the uncertainty of dedicated VLBI Cali- 
bration Surveys. 

We thank an anonymous referee for constructive sugges- 
tions. The National Radio Astronomy Observatory is a facil- 
ity of the National Science Foundation operated under coop- 
erative agreement by Associated Universities, Inc. 

Facilities: VLBA (project code BT085). 
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